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Profiling, bioinformatic, and functional data on the 
developing olfactory/GnRH system reveal cellular and 
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potentially relevant for the Kallmann syndrome 
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During embryonic development, immature neurons in the olfactory epithelium (OE) extend 
axons through the nasal mesenchyme, to contact projection neurons in the olfactory bulb. 
Axon navigation is accompanied by migration of the GnRH+ neurons, which enter the ante- 
rior forebrain and home in the septo-hypothalamic area. This process can be interrupted at 
various points and lead to the onset of the Kallmann syndrome (KS), a disorder characterized 
by anosmia and central hypogonadotropic hypogonadism. Several genes has been identi- 
fied in human and mice that cause KS or a KS-like phenotype. In mice a set of transcription 
factors appears to be required for olfactory connectivity and GnRH neuron migration; thus 
we explored the transcriptional network underlying this developmental process by pro- 
filing the OE and the adjacent mesenchyme at three embryonic ages. We also profiled 
the OE from embryos null for Dlx5, a homeogene that causes a KS-like phenotype when 
deleted. We identified 20 interesting genes belonging to the following categories: (1 ) trans- 
membrane adhesion/receptor, (2) axon-glia interaction, (3) scaffold/adapter for signaling, (4) 
synaptic proteins. We tested some of them in zebrafish embryos: the depletion of five (of 
six) Dlx5 targets affected axonal extension and targeting, while three (of three) affected 
GnRH neuron position and neurite organization. Thus, we confirmed the importance of 
cell-cell and cell-matrix interactions and identified new molecules needed for olfactory 
connection and GnRH neuron migration. Using available and newly generated data, we 
predicted/prioritized putative KS-disease genes, by building conserved co-expression net- 
works with all known disease genes in human and mouse. The results show the overall 
validity of approaches based on high-throughput data and predictive bioinformatics to iden- 
tify genes potentially relevant for the molecular pathogenesis of KS. A number of candidate 
will be discussed, that should be tested in future mutation screens. 
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INTRODUCTION 

Central Hypogonadic Hypogonadism (CHH), is a heterogeneous 
genetic disorders characterized by absent or incomplete puberty, 
due to low circulating gonadotropins and sex steroids. Its mode of 
inheritance can be X-linked, autosomal dominant, or autosomal 
recessive, although unrelated sporadic cases occur more frequently 
(1). The disease is often associated with anosmia/hyposmia, in this 
case it is known as Kallmann Syndrome [KS, on-line Mendelian 
inheritance in man (OMIM) 308700], or with a normal sense of 
smell (normosmic CHH, or nCHH). These conditions are variably 
associated with non-reproductive phenotypes such as unilateral 
renal agenesis, skeletal abnormalities, midline malformations, or 
hearing loss. Neurological symptoms (including synkinesia of the 



hands, sensorineural deafness, eye-movement abnormalities, cere- 
bellar ataxia, and gaze-evoked horizontal nystagmus) may also 
occur depending on the specific mode of inheritance (2). 

Mutations affecting a large number of unrelated genes 
have been linked to the onset of KS/nCHH, currently includ- 
ing Anosminl (KALI), Fibroblast Growth Factor Receptor-1 
(FGFR1), Fibroblast Growth Factor 8 (FGF8), GnRH receptor 
(GNRH-R), Nasal Embryonic LHRH Factor (NELF), Kisspeptin 
(KISS1); Kisspeptin Receptor (KISS-R)/G-protein-Coupled Recep- 
tor 54 (GPR54), Prokineticin-2 (PROK-2), Prokineticin Receptor- 
2 (PROKR2), Chromodomain Helicase DNA-binding Protein 7 
(CHD7), Neurokinin-B (TAC3), Neurokinin-B Receptor (TAC3R), 
Heparan Sulfate 6-O-SulphoTransferase 1 (HS6ST1 ), SOX10, 
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Semaphorin-3A (SEMA3A), and five novel genes, members of the 
"FGF8-synexpressome" (1-8). In addition, several mouse models 
of targeted gene disruption have been shown to exhibit a KS-like 
phenotype (6, 9-18). 

Despite the number of genes mutated in KS/nCHH, the major- 
ity of patients (>60%) do not harbor mutations in known disease 
genes, thus it is expected that many additional disease loci remain 
to be identified. In addition, the mutations found in KS/nCHH 
patients, once thought to act alone, are now recognized as cooper- 
ating mutations, and in fact in some cases a bi- genie or oligo- genie 
origin of these disease has been reported, with specific geno- 
type/phenotype correlations (19-22). These findings open ques- 
tions on the actual prevalence of single and combined mutations, 
the functional cooperation between them, and the possibility to 
use these information for accurate prognostic evaluations. 

Kallmann syndrome is righteously considered a developmental 
disease. During embryonic development the GnRH neurons orig- 
inate in the primitive olfactory area, migrate along the extending 
axons from the olfactory epithelium (OE) and the vomero-nasal 
organ (VNO), reach the anterior-basal forebrain and home in 
the septo-hypothalamic region of the adult brain (23-28). The 
association of the olfactory axons with the immature GnRH neu- 
rons, hence their ability to migrate and reach the hypothalamus, 
is an ancient and highly conserved developmental process, justi- 
fied by the fact that it is essential for puberty and reproduction, 
in addition to neuro-modulatory functions (27, 29, 30). Not sur- 
prisingly this process is governed by a large set of molecular cues. 
Several studies have identified specific signaling molecules and 
their cognate receptors, as well as adhesion molecules, axon-glia 
and axon-matrix molecules play a role in guiding the axons to 
the correct position and consent the penetration of the basement 
membrane and the brain parenchyma (31-40). For instance, the 
semaphorin co-receptor Nrpl is expressed by extending axon and 
GnRH neurons, and mediates the guiding functions of Sema3a, 
expressed in the nasal mesenchyme ( 10, 1 1 ). FGF8 has been shown 
to act as survival factor for olfactory and migrating GnRH neu- 
rons, which express its receptor FGFR1 (41-43), and both genes 
are mutated in a subset of KS/nCHH patients. To further com- 
plicate the picture, a cell population on the surface of the OB 
interacts with incoming axons, GnRH neurons and the CNS, and 
provide key signals for basement membrane fenestration, hence 
axon connectivity (44). 

Due to the close relationship of olfactory axon elonga- 
tion/connectivity and GnRH migration that occurs during embry- 
onic development, the GnRH neuronal migration is strictly depen- 
dent on the integrity and connectivity of the olfactory pathway 
(30). A premature termination or mislocalization of olfactory 
axons results in impaired odor perception and GnRH homing. 
Thus, defects in olfactory development and/or GnRH neuron 
migration are considered the main primary cause of KS. The 
genetic findings summarized above have revealed much about 
the abnormalities that can befall both the development of the 
olfactory sensory system and GnRH neuron ontogenesis, includ- 
ing their differentiation, migration, maturation, circuit formation, 
and senescence. 

Experimentally, animal models with altered olfactory and 
GnRH development are becoming available, including mouse, 



Zebrafish, and Medaka. The zebrafish embryo is ideal for develop- 
mental genetic studies, and the depletion of anosmin-la leads to 
altered olfactory development and a KS-like phenotype (45, 46). 
In mice, several mutant strains display a phenotype that closely 
resemble KS/nCHH, including mouse mutant for Dlx5 (14, 16, 
47), Emx2 (18), Klf7 (13), FezFl (17, 48), Sixl (12), Prok2 and its 
receptor Prokr2 (6, 15), Lhx2 (9), Ebf2 (49), Nrpl and Sema3a (10, 
11). Notably, 7 of these (Dlx5 y Emx2 y Klf7 y Sixl, FezFl, Ebf2 y and 
Lhx2) code for unrelated transcription factors, thus it can be pos- 
tulated the existence of transcription regulatory networks, yet to 
be uncovered, that sustain olfactory development and connectivity, 
consent migration of the GnRH neurons and may contribute to the 
onset of KS/nCHH when altered. Furthermore, it is increasingly 
recognized that biological processes are governed and regulated 
by regulatory modules and networks of molecular interaction, 
not limited to protein-coding genes, rather then simplistically by 
individual genes. 

To advance in our knowledge on the molecular regulation of 
axon extension/connectivity and GnRH neuron migration, in the 
present study we adopted a strategy based on the generation of 
transcriptome-wide profile data, combined with bioinformatic 
analyses and meta-analyses. In addition to the normal olfactory 
tissue we have also included one of the mouse models of KS, i.e., 
the Dlx5 null (14, 16, 47). We then used transgenic Zebrafishes 
to image the olfactory axons and the GnRH neurons, and use 
these to establish the function of Dlx5 targets for olfactory axon 
extension/contact and on GnRH neuron migration and neurite 
extension. The results confirm a role for Dlx5 and FGFR1, and 
indicate Lrrnl and Lingo2 as novel players for olfactory axon orga- 
nization and for GnRH neuron migration. Finally, we applied 
a gene prediction algorithm based on conserved co-expression 
networks, on all known human and mouse KS-causing genes. 
We predict a set of best candidates for causing, con-causing, or 
modifying the KS/nCHH phenotype. 

MATERIALS AND METHODS 
MICE NULL FOR Dlx5 

Mice with targeted disruption of Dlx5 have been previously 
reported (50). The null allele, denominated Dlx5 lacZ , allows 
for detection of the D/x5-expressing cells by staining for |3- 
galactosidase (|3-gal) expression. The olfactory phenotype has been 
previously characterized (14, 16, 47). To obtain the WT sam- 
ples, only WT males and females were crossed. To obtain Dlx5 
mutant samples, Dlx5 +I ~ (heterozygous) males and females were 
crossed; the progeny showed the expected Mendelian ratios of 
genotypes +/+ , Dlx5 +I ~ and Dlx5~ l ~ . Pregnant females were sac- 
rificed at the chosen embryonic age by cervical dislocation. The 
day of the vaginal plug was considered E0.5. All animal proce- 
dures were approved by the Ethical Committee of the University 
of Torino, and by the Italian Ministry of Health. 

TISSUE COLLECTION FROM MOUSE EMBRYOS 

Embryos were collected clean of extra- embryonic tissues (used 
for genotyping) by manual dissection, transferred in RNAse-free 
PBS, and further dissected to separate the head. This was then 
included in 3% low-melting agarose in PBS, let harden and sec- 
tioned by vibratome (250 |xm). Sections were manually dissected 



Frontiers in Endocrinology | Experimental Endocrinology 



December 2013 | Volume 4 | Article 203 | 2 



Garaffo et al. 



Olfactory and GnRH development 



in cold PBS, with fine pins, to collect the OE or the VNO epithe- 
lia, or alternatively to collect the adjacent mesenchyme (Figure SI 
in Supplementary Material). The excised tissues were individu- 
ally collected in RNA-later (Ambion) and stored at — 20°C until 
extraction. Following genotyping, samples of the same genotype 
were pooled. For the Dlx5 mutant tissues, the entire epithelial lin- 
ing of the nasal cavity was collected, since it was not possible to 
discriminate the OE vs. the respiratory epithelium. 

RNA EXTRACTION, LABELING, AND HYBRIDIZATION ON MOUSE 
EXON-SPECIFIC ARRAYS 

At least 15 embryos were used for each developmental age, the col- 
lected tissues were pooled in three independent biological samples, 
used to extract total RNA with the Trizol (Invitrogen). 

After extraction, RNA samples were quantified using a 
NanoDrop spectrophotometer (NanoDrop Technologies), the 
integrity of RNA molecules was assessed by capillary electrophore- 
sis on a Agilent Bioanalyzer (Agilent), and found to have a RIN 
(RNA Integrity Number) value >5. One microgram of each total 
RNA sample (in triplicate) was processed using the Affymetrix 
platform's instruments, following the GeneChip Whole Transcript 
Sense Target Labeling procedure, according to instructions. Ribo- 
somal RNA was depleted using the RiboMinus kit (Invitrogen), 
cDNA was synthesized with random primers coupled with the 
T7 Promoter sequence, using Superscript II for first-strand syn- 
thesis, and DNA Polymerase I for second-strand synthesis. The 
cDNA was and used as template for IVT amplification, using 
T7 polymerase. The amplificated products were used to syn- 
thesize single-stranded cDNAs, with the incorporation of dUTP, 
the products were fragmented by uracil- DNA- glycosylase (UDG) 
and apurinic/apyrimidinic endonuclease-1 (APE 1) treatment. 
Finally, 5.5 |xg of fragmented cDNA samples were biotinylated 
with terminal deoxynucleotidyl transferase and used to hybridize 
on GeneChip® Exon 1.0 ST Arrays (Affymetrix, Santa Clara, 
USA). The Chips were washed and stained with Streptavidin- 
phycoerithrin in the GeneChip Fluidic Station 450 and scanned 
with Affymetrix GeneChip® Scanner 3000 7G. 

ANALYSIS OF MICR0ARRAY DATA 

Quality control was performed using the Affymetrix Expression 
Console software 1 . All the experiments exhibited optimal qual- 
ity controls and correctly clustered in the right sample groups; 
they were thus all included in the analysis. Normalization and 
probeset summarization steps were performed with RMA, within 
the OneChannelGUI package (51) included in Bioconductor (52), 
separately for each pairwise comparison including the six rele- 
vant arrays (three biological replicates per condition). Differen- 
tially expressed genes (DEG) for each pairwise comparison were 
obtained with Rank Products (53), adopting a 0.05 false discovery 
rate (adj. p-value <0.05). 

SOFTWARES AND DATABASES 

For preliminary Gene Ontology (G.O.) analyses we used DAVID 2 
and KEGG 3 . For improved categorization and visualization, we 



1 www.affymetrix.com 

2 http://david.abcc.ncifcrf.gov/ 

3 http://www.genome.jp/kegg/pathway.html 



used ClueGO (54). For the time course analysis we used default 
parameters. For the analysis of down-regulated DEGs in the 
Dlx5~ f ~ samples we relaxed the analysis by using a cutoff of 
0.001 on nominal enrichment p- value. For embryonic expression 
of RefSeq genes we used the two on-line in situ hybridization 
databases GenePaint 4 and Eurexpress 5 . For the position weight 
matrix (PWM) we used the JASPAR database. Tissue-specific con- 
served co-expression networks were obtained with the TS-CoExp 
Browser 6 (55). 

We also used the following web resources: Ensembl Genome 
Browser 7 , UCSC Genome Browser 8 , RefSeq 9 , Mouse Genome 
Informatix 10 ,OMIM n . 

GENOME-WIDE PREDICTION OF Dlx BINDING SITES AND PUTATIVE 
TARGET GENES 

With the PWM of Dlx5 provided by JASPAR under acces- 
sion PH0024.1 (56) Dlx5 sites were predicted by standard log- 
likelihood ratios, using as null model the nucleotide frequen- 
cies computed over the whole intergenic fraction of the mouse 
genome. We considered only those sites scoring 50% of the 
maximum possible score or better. We selected sites that are 
conserved in at least two (of eight vertebrate species). A site 
is defined as conserved with species S if it lies in a region 
of the mouse genome which is aligned with a region of the 
S genome and the aligned sequence in/S/is a site according 
to the same definition used for mouse sites. A ranked list of 
putative Dlx5 targets was obtained from the identified sites as 
described (57). 

CONSERVED CO-EXPRESSION NETWORK, AND 
PREDICTION/PRIORITIZATION OF PUTATIVE DISEASE GENES 

Tissue-specific conserved co-expression networks were obtained 
with the TS-CoExp Browser (see footnote text 6) (55, 58), based 
on 5188 human and 2310 mouse manually annotated microar- 
ray experiments. For disease prediction/prioritization we used 
a tool within the TS-CoExp Browser and the same approach 
based on conserved co-expression networks, but instead of 
using genes causing similar phenotypes, we used KS-disease 
genes as "reference" genes. These genes were selected based 
on documented mutations in KS patients (for human) or well 
described olfactory/GnRH embryonic phenotype recapitulating 
KS (mouse). 

VALIDATION OF ARRAY DATA BY REAL-TIME qPCR 

Tissue samples corresponding to WT and Dlx5~ l ~ OE were col- 
lected from embryos at the age El 2.5, transferred in RNA-later 
in individual tubes and stored at — 20°C. The genotype was 
determined on extra- embryonic tissues. Samples were pooled 
according to the genotype, collected in Trizol (Invitrogen), and 



4 www. genepaint. org 

5 www.eurexpress.org 

6 http ://www.mbcunito.it/cbu/ts - coexp 

7 http://www.ensembl.org/index.html 

8 http://genome.ucsc.edu 

9 http://www.ncbi.nlm.nih.gov/RefSeq 
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used to extract total RNA according to the instructions. For 
Real-Time qPCR, 250 ng of total RNA was reverse-transcribed at 
42°C for 50 min in the presence of 500 ng/|xl random hexamers, 
10 mM of each dNTPs, RNasin and Improm Reverse Transcrip- 
tase (Promega). Relative cDNA abundance was determined using 
the AB7900 System and the GoTaq qPCR Master Mix (Promega). 
Specific cDNAs were amplified using primers and probes designed 
according the Universal Probe Library system (UPS, Roche). 
Experiments were repeated at least twice on independent sam- 
ples, every point was done in triplicate, results were normal- 
ized to the level of TATA-binding protein (TBP) and GAPDH 
mRNAs. Data analysis was performed with ABI software, ver- 
sion 2.1 (Applied Biosystems) using the comparative Cq method, 
calculated with the formula of the DDCq. For each primer-pair, 
the melting curves of the amplified products revealed a single 
peak. Primer sequences are provided (Table SI in Supplementary 
Material). 

ZEBRAFISH STRAINS AND GENE KNOCK-DOWN IN EMBRYOS 

The following two strains were used for visualization of the olfac- 
tory axons: OMP 2k :gap-CFP rw0M and TRPC2 4 - 5k :gap-Venus rw037 
(59-61), and were obtained from Drs. Nobuhiko Miyasaka and 
Yoshihiro Yoshihara (RIKEN Brain Science Inst., Japan). The fish 
strain GnRH3:GFP (62-64) was obtained from Dr. Y. Zohar (Uni- 
versity of Maryland Biotechnology Institute, Baltimore, USA) and 
Dr. Y. Gothilf (Life Sciences, Tel- Aviv University, Israel). Adult 
fishes were maintained, bred and genotyped according to standard 
procedures, kept under a 14 h-light and 10 h-dark photoperiod at 
28°C. Allelic transmission followed the expected Mendelian ratios. 
Following fertilization, 1-cell zygotes were collected and main- 
tained in the presence of 0.003% l-phenyl-2-thiourea (PTU) to 
prevent formation of melanin. 

To down-modulate specific genes, we injected antisense mor- 
pholino oligos (MO) into zebrafish oocytes (65, 66). MO were 
designed either to block splicing at a specific exon-intron junc- 
tion (GeneTool oligo design), and consequently lead to present of 
aberrant transcripts and frame-shifted translation, or to anneal 
to the ATG start codon and inhibit translation initiation. For 
z-dlx5a we combined two MOs: one annealing with the exonl- 
intronl splice junction and leading to a premature Stop codon 
upstream of the homeodomain; the other annealing with the Start 
codon. Sequences and properties of all the MO are in Table S2 in 
Supplementary Material. Zygotes were collected at one-cell stage 
and injected under stereological examination with 4 ng of MO, 
in presence of Phenol Red for subsequent selection. From 48 to 
72 h post fertilization (hpf) embryos were fixed with 4% PFA 
at 4°C ON, washed in PBS, and embedded in 4% low-melting 
agarose, 0.1% Tween-20. The apical portion of the head was man- 
ually dissected from the rest of embryo. Confocal microscopy 
analysis was performed using a Leica TCS SP5 (Leica Microsys- 
tems). The OMP:CFP+ and the Trpc2:Venus+ (YFP+) axons 
were viewed in a frontal plane, while the GnRH3:GFP+ neurons 
were viewed in a ventral plane. Images were acquired as Z-stacks 
of 1 |xm thick optical sections. Digital micrographs images were 
contrast balanced and color matched using Photoshop7 (Adobe), 
cropped, rotated, and assembled into figures with QuarkXpress 
(Pantone). 



RESULTS 

GENES DIFFERENTIALLY EXPRESSED DURING OLFACTORY 
DEVELOPMENT 

We set forth to generate expression profiles of the OE at key stages 
of its development, comprising the time of axonal connection. 
We selected three developmental stages, i.e., the Olfactory Pla- 
code (OPL) at El 1.5, the OE at E12.5, and either the OE or the 
VNO at E14.5. Mouse Affymetrix GeneChip® Exon 1.0 ST Arrays 
were used to analyze the gene expression profiles of the develop- 
ing olfactory (neuro) epithelium (OE). Comparing the OE E12 vs. 
the OE at Ell, with adj. p-value <0.05 and fold-change <-0.9 
or >0.9, we found 29 up-regulated and 62 down -regulated genes. 
Comparing the OE at E14 vs. the OPL Ell we found 358 up- 
regulated and 17 down -regulated genes. Comparing the VNO E14 
vs. the OPL El 1 we found 459 up-regulated and 2 1 down-regulated 
genes. 

A fraction of the DEGs might derive from mesenchymal cells 
present in the epithelial samples; as a matter of fact, epithelial 
cells do not easily detach from the basement membrane and 
mesenchymal cells inevitably tend to remain attached. A survey 
of the embryonic expression territory of the modulated genes 
using the on-line expression databases and www.genepaint.org 
and www.eurexpress.org showed that about 10% of the DEGs was 
indeed expressed in the nasal mesenchyme adjacent to the OE, 
and not in the OE or VNO proper. Thus, we decided to estimate 
the extent of mesenchymal contamination in the OPL, OE, and 
VNO samples, by collecting pure mesenchymal tissue adjacent to 
the OPL, OE, and VNO, at the same embryonic ages, and use the 
RNA extracted from these to quantitatively determine the mRNA 
abundance of "epithelial only" (Foxjl, Fmo2, and Ehf) and "mes- 
enchymal only" (Sp7 and Lectl) genes, by Real-Time qPCR. In 
the same experiment we compared the samples of the OE (mixed 
epithelium and mesenchyme) with "pure mesenchyme" samples at 
the same embryonic age. The results indicate that the abundance of 
a mesenchymal RNA in the OE samples is roughly 15% that of the 
pure mesenchyme samples, thus we assumed that the contribution 
of MES in the EPI samples is 15% (Figure S2 in Supplementary 
Material). 

At the same time, using the GeneChip® Exon 1.0 ST Arrays 
and the same hybridization procedure and statistical analyses 
used before, we generated profiling data from the MES sam- 
ples collected from wild-type embryos at El 1.5, E12.5, and 
E14.5. At the age El 4.5 the samples were collected adjacent to 
the OE or adjacent to the VNO, according to their anatomical 
position, and maintained separated. This effort was undertaken 
to: (1) explore the global changes of expression that under- 
lie interaction between the OE and the MES, (2) carry out 
a subtraction step on the raw EPI data, to generate cleaner 
OE data. 

By comparing the MES samples at E12 vs. Ell we found 118 
up-regulated and 17 down -regulated genes; comparing the sam- 
ples OE at E14 vs. OPL Ell we detected 284 up-regulated and 
41 down-regulated genes, while comparing the VNO at E14 vs. 
OPL Ell we detected 293 up-regulated and 35 down-regulated 
genes (the non-annotated probes are not included). Then we sub- 
tracted the estimated expression of MES genes (15%) from the raw 
expression data, applying this general formula to all genes present 
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and expressed: 

EHs)=Ei(g)-FxM(g) 

where is the expression of gene g in the i-th replicate of 

the EPI dataset and M(g) is the expression in the MES dataset, 
averaged over all replicates. F is the estimated mesenchymal frac- 
tion, equal to 0.15. Choosing F to be equal to 0.1 or 0.2 did 
not significantly alter the results. With this calculation we created 
a subtracted and corrected dataset with expression values more 
indicative of the sole EPI expression. Comparing the corrected 
EPI samples at E12 vs. Ell we found 9 up-regulated genes and 
57 down-regulated; comparing the samples OE at E14 vs. OPL 
Ell we detected 250 up-regulated and 19 down -regulated genes, 
while comparing the VNO E14 vs. OPL Ell we detected 347 up- 
regulated genes and 14 down-regulated (the non-annotated probes 
and the OR genes are not counted). After the subtraction, a num- 
ber of genes reached a "no expression" level. We assume that this is 
due mainly to the fact that their differential expression was relative 
to the MES. We examined how many of the genes that disappeared 
from the raw list are up-regulated in the MES samples, and detected 
highly significant enrichments (p < 4e-12). 

The corrected lists of EPI DEGs up- and down-modulated in 
the OE E14 vs. OPL Ell are reported in Tables S3 and S4 in Sup- 
plementary Material, respectively, while the lists of DEGs up- and 
down-modulated in the VNO E14 vs. OP Ell are in Tables S5 



and S6 in Supplementary Material, respectively. In the OE we find 
genes expected to be associated or to play a role in neuronal dif- 
ferentiation and/or olfactory development, such as NeuroD, OMP, 
Peripheric NCAM2, Claudius, Keratins, and Lhx2 (a gene caus- 
ing a KS-like phenotype in the mouse) (9). In addition we find 
a set of olfactory receptor (OR) genes, as expected (Table S7 in 
Supplementary Material). 

Next we carried out functional categorization analyses on the 
genes up-regulated in the OE, to identify enriched functional cate- 
gories, using the Gene Ontology-based ClueGo tool (54, 67). Since 
this analysis could be biased by the OR genes, which are numerous 
(about 1000 in the mouse genome) and belong to a single cate- 
gory, we masked the OR genes. The results are shown in Figure 1 A. 
From the comparisons OE 14 vs. OPL 11 we detect: regulation of 
epithelial cell proliferation, regulation of cell migration, regula- 
tion of extracellular matrix organization, and various categories 
of response to signals. 

In the VNO, we find several genes expected to be associ- 
ated or to play a role in VNO development, such as NeuroD, 
OMP, Lhx2, Peripheric Claudius, Keratius, EphA3, Neuropiliul, 
Lamiuiu$3, Lhx2 (a Kallmann gene in the mouse), and Dcx. In 
addition we find several OR genes, as expected (Table S8 in Sup- 
plementary Material). We carried out functional categorization 
on the genes up-regulated in the VNO, after masking the differ- 
entially expressed OR genes, and detected the over-represented 
classes shown in Figure IB. Focusing on the comparisons E14 vs. 




FIGURE 1 | Continued 
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FIGURE 1 | Functional categorization of genes differentially 
expressed during normal OE and VNO development. Significantly 
over-represented functional categories of DEGs detected in the 
(corrected) EPI samples, comparing OE 14.5 vs. OPL 11.5 (A) or 



El 1, we detect: regulation of epithelial cell proliferation, regulation 
of cell migration, regulation of cell adhesion, gland and epithelium 
morphogenesis, cartilage development, bone development, extra- 
cellular matrix organization, and various categories of response to 
signals. 

DEGs IN OLFACTORY-ASSOCIATED MES, DURING DEVELOPMENT 

We then compared the profiles of the MES samples across the 
developmental ages E11.5-E12.5-E14.5. The full lists of up- and 



comparing the VNO 14.5 vs. OPL 11.5 (B) were organized them in 
networks using ClueGO. Only the results of the up-regulated DEGs are 
shown. Circles represent over-represented categories, lines indicate 
related categories. 



down-modulated DEGs relative to the OE are provided in 
Tables S9 and S10 in Supplementary Material, respectively, while 
the full lists of up- and down- regulated DEGs relative to the 
VNO are provided in Tables SI 1 and S12 in Supplementary Mate- 
rial, respectively. We recognized genes playing a role in cell-cell 
communication, signaling, matrix remodeling, etc. such as Inte- 
grins, Contactins, Matrillins, Tenascin, Collagens, MMPs, Adams, 
Lectin Galactose Binding 9, Elastin, FGF7, FGF12, Sfrp2, Sfrp4, 
Sema3D, Sema3C, Nrpl, Wnt2, Bmp5, Follistatin. We also found 
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some neuronal genes, likely due to a minimal presence of olfactory 
neuron in the MES sample and to the presence of migratory GnRH 
neurons in the E14 sample, minimal in the Ell sample. Functional 
categorization on these DEGs detected an enrichment of the fol- 
lowing categories: extracellular matrix organization, cell- substrate 
adhesion, cartilage and bone development, organ morphogenesis, 
response to signals, and some neuronal categories (Figures 2A,B). 

PROFILING OF THE DIX5" 7 - VS. WILD-TYPE OE 

The Dlx5~ l ~ mutant mice represent a fully penetrant model of 
KS (14, 16, 47). Triplicates of the OE and VNO tissues were col- 
lected from WT and D/x5 _/ ~ embryos at the age E12.5, total 
RNA was extracted and hybridized on the GeneChip® Exon 1.0 
ST Arrays. Using the indicated statistical parameters (see Materi- 
als and Methods) we detected 121 down- and 25 up-regulated 
genes in the Dk5 _/_ OE vs. the WT, not counting the non- 
annotated probes and the OR genes (Figure 3A; Table S13 in 
Supplementary Material). Again, the OR genes were removed 



(provided in Table S14 in Supplementary Material) prior to con- 
ducting functional categorization analysis by G.O. We detected: 
intermediate filament/cytoskeletal organization, endocrine sys- 
tem development, forebrain development, cell-cell signaling, and 
epithelial cell differentiation (Figure 3B). 

We carried out a technical validation of the microarray results, 
by selecting 12 down- and 4 up-regulated DEGs and quantifying 
their expression on independently collected samples, by Real- 
Time qPCR. Of these DEGs, 1 1 down- and 2 up-regulated were 
confirmed (Figure 3C). 

Next we verified whether the identified DEGs are expressed in 
the OE, in the adjacent mesenchyme, or in both, by consulting the 
on-line expression databases GenePaint and Eurexpress. We classi- 
fied DEGs as either not expressed (— ), expressed in the neurepithe- 
lium (NEp), expressed in the mesenchyme (Mes), expressed in the 
respiratory epithelium (Res) or ubiquitously expressed (Ub). We 
assumed that the OR were all expressed in NEp, and in any case they 
were excluded. Among the 72 down-modulated DEGs considered, 
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FIGURE 2 | Functional categorization of genes differentially expressed in 
the mesenchyme adjacent the OE and VNO, during development. 

Significantly over-represented functional categories of the DEGs detected in 



the mesenchymal tissue associated with the OE 14.5 vs. OPL 11.5 (A) or 
associated with the VNO 14.5 vs. OPL 11.5 (B) were organized in networks 
using ClueGO. Only the results of the up-regulated DEGs are shown. 



50 (69%) have a NEp expression, 14 (19%) are not expressed in the 
OE, 4 (5%) have a Mes expression, 3 (4%) have a Res expression, 
and 2 (3%) are ubiquitously expressed. Thus, conclude that the 
majority of down -modulated genes are expressed in the OE. 

Then we examined whether the DEGs were differentially 
expressed also in other mouse tissues upon loss of Dlx5, specifically 
the inner ear and the pharyngeal arches (68, 69). No common DEG 
was found, indicating that Dlx5 targets are strongly tissue-specific. 
Next we examined whether the olfactory DEGs we detected were 
also differentially expressed in other published mouse models of 
KS, i.e., the Klf7 and the Emx2 mutants (70, 71). Three genes were 
found commonly differentially expressed the three models, namely 
stathmin-like 3, synaptotagmin i, and calmegin, all expressed in 
the embryonic OE. Fifty genes were in common between Dlx5 
and Emx2 datasets, seven were in common between Emx2 and 
Klf7 y and one was in common between Dlx5 and Klf7. However it 



should be noted that the profiles of the Emx2 and Klf7 mutants 
were obtained from the OB and not the OE. 

The DEGs up -regulated in the absence of Dlx5 are enriched 
in generic terms: biosynthesis, metabolic processes, morphogen- 
esis. Of the 27 DEGs considered, 14 (52%) are expressed in the 
Mes, 6 (22%) are not expressed in the OE/VNO, 3 (11%) have a 
Res expression, 2 (7.5%) are expressed in the OE, and 2 (7.5%) 
are ubiquitously expressed. Thus, we conclude that most of the 
up-regulated DEGs are not expressed in the OE. Since the Dlx pro- 
teins are generally considered transcription activators (72, 73), the 
interest in these DEGs is low and they were not further considered. 

GENOME-WIDE PREDICTION OF Dlx5 BINDING SITES AND 
TRANSCRIPTIONAL TARGETS 

Using the consensus PWM for Dlx5 (56) (Figure S3 in Supple- 
mentary Material) we screened conserved regions of the vertebrate 
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B 



Term P value 

GO:0045104~intermediate filament cytoskeleton 3 E-07 

GO:0045103~intermediate filament-based process 6 E-07 

GO:0045109~intermediate filament organization 8.1 E-05 

GO:0007267~cell-cell signaling 7.7 E-04 

GO:0019226~transmission of nerve impulse 8.4 E-04 

GO:002 1 983~pituitary gland development 9.6 E-04 

GO:0035270~endocrine system development 0.00142 

GO:00070 1 0~cytoskeleton organization 0.00 1 63 

GO:002 1 536~diencephalon development 0.002 1 3 

GO:0030900~forebrain development 0.00549 

GO:0007268~synaptic transmission 0.007 1 6 

GO:0030855~epithelial cell differentiation 0.00968 

GO:0048732~gland development 0.0 1 083 

GO:003 1 424~keratinization 0.0 1 528 




FIGURE 3 | Genes differentially expressed in the D\x5~'~ OE vs. 
wild-type OE. (A) Heatmap and cluster analyses on the MicroArray 
data, after elimination of non-annotated genes and of the OR 
genes. The triplicate Dlx5~'~ (KO) and wild-type (WT) samples are 
indicated at the bottom. (B) Gene Ontology to detect 



over-represented functional categories. (C) Validation by Real-Time 
qPCR on selected down-regulated genes, done on independently 
collected samples. The abundance in the wild-type tissue is set = 1. 
For each gene, the fold-change determined by MicroArray 
hybridization is also reported. 
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genome and detected putative Dlx5 binding sites. We attributed to 
each site a score that reflects the number of species in which the site 
is conserved. We then associated the sites to an associated Refseq 
transcripts and found 3,426 RefSeq targets, corresponding to 2,683 
unique Entrez-IDs [see Materials and Methods, and Ref. (57)]. 
The top scoring RefSeq are reported in Table S15 in Supplemen- 
tary Material. We then categorized the predicted Dlx5 targets by 
ClueGO (54), and detected an enrichment categories such as neu- 
ronal differentiation, brain development, etc. as expected [there is 
ample literature on this; see Ref. (73)]. 

We then intersected the best predicted Dlx5 targets (having 
at least one binding site conserved in at least three mammalian 
species, and located <10kb from the TSS) with the list of DEGs 
obtained comparing Dlx5~'~ vs. WT OE, and we found that 16% 
of the down-regulated DEGs (19/121; p = 0.0003) were indeed 
predicted target of Dlx5, while 40% of up-regulated DEGs (9/21; 
p = 0.00019) were predicted targets. In both cases statistical signif- 
icance was reached. This suggests that the prediction algorithm we 
have used is sensitive and sound. To restrict the number of candi- 
date genes we intersected the profile datasets with: (a) embryonic 
expression databases, (b) conserved co- regulations, (c) predicted 
Dlx5 sites and target RefSeq, (d) data from published literature, 
in order to assign a score value to each DEG (Tables 1A-E). The 
expression of these putative Dlx5 targets in the embryonic OE 
and nasal region, by in situ hybridization (see footnote text 4) 
is reported in Figure S5 in Supplementary Material. Some of the 
most functionally relevant genes are briefly described below: 

Lrrnl codes for a transmembrane protein related to Drosophila 
TRN/CAPS proteins, known play a role in neuromuscular tar- 
get recognition, and to mediate interactions between incoming 
axons and the targets, possibly via homophilic adhesion. Lrrnl is 
expressed in the mouse embryonic OE. 

Lingo2 (also known as Lrrn6c) codes for a transmembrane pro- 
tein, expressed in the OE and in the ventricular region of the 
embryonic forebrain. Lingo proteins interact with the NOGO 
receptor and are able to modulate the NOGO pathway (74), 
however their precise functions are poorly known. 

Lgil codes for a leucine-rich repeat secreted molecule of the 
SLIT family, involved in growth of neuronal processes on myelin 
substrates (75, 76). 

St8siaVI is expressed by olfactory neurons and might be 
implicated in polysialylation the N-CAM to confer anti-adhesive 
properties to neuronal surfaces (77-79). 

Homer2 codes for a protein present at post-synaptic density, 
likely to be involved in receptor clustering and trafficking, as well 
as calcium homeostasis (80). Recently, a role of Homer2 in tun- 
ing the activity of G protein-coupled receptors (such as ORs) has 
been reported (81). Homer2 is expressed in the OE of the mouse 
embryo, however its function is unknown. 

TESTING Dlx5, Dlx5 TARGETS, AND KS GENES IN ZEBRAFISH STRAINS: 
THE OLFACTORY AXONS 

The development of olfactory system is well conserved during 
vertebrate evolution (27, 35, 82) and consists of two indepen- 
dent components: the main OE for detecting chemical com- 
pounds (odorants) and the VNO-accessory system for detecting 



pheromones. Fishes and primates lack a VNO organ and present 
only one olfactory organ, the OE (83). Within the OE of the fish, 
all ORNs project their axons to the OB - at different region in a 
mutually exclusive manner (60) - but display distinct properties 
with respect to their morphology, relative position in the OE, and 
molecular expression. The ciliated OSNs with long dendrites are 
situated in the deep layer of the OE, whereas microvillous ORNs 
with short dendrites are located in the superficial layer. The cili- 
ated and microvillous ORNs are reported to express OR-type and 
V2R-type receptors, respectively (84, 85). 

We opted to use Danio rerio (zebrafish) as a model to func- 
tionally examine in vivo the identified DEGs for their role in 
olfactory/GnRH development. We used two transgenic zebrafish 
strains expressing distinct fluorescent proteins in the fish olfac- 
tory neurons (59-61). In one strain the CFP reporter is expressed 
under the control of OMP promoter, which marks the majority 
of basal-layer ORN, projecting their axons to the dorsal OB. In 
the other strain, the Venus (YFP) reporter is expressed under the 
control of the Trpc2 promoter, which marks a sub-population of 
apical-layer ORN, projecting to the ventro -lateral OB (scheme in 
Figure S4 in Supplementary Material) (60). The CFP+ and the 
Venus+ (YFP+) neurons are thought to correspond, respectively, 
to the OE and VNO receptors of the mammalian system. Since 
the reporter fluorescent proteins are efficiently translocated in the 
ORN axons, these two strains visualize the peripheral olfactory 
pathway. 

We tested z-fgfrla, the fish ortholog of mammalian FGFR1, 
to establish whether its depletion recapitulates the hallmarks of 
KS. Notably, mice hypomorphic for FGF8 expression show dis- 
tinctive signs of a KS phenotype, i.e., impaired migration of 
GnRH+ neurons and defects in olfactory development (41, 42). 
We injected z-fgfrla MOs in 1-cell embryos of the OMP:CFP and 
the Trpc2:Venus strains, and 72hpf we examined the number of 
fluorescent embryos, the organization of the OPL, the fascicula- 
tion, extension and glomeruli formation. In 61% (32/52) of the 
embryos we observed an altered morphogenesis of the OPL and 
an abnormal distribution of the CFP+ and the Venus+ neurons 
within the OLP (Figure 4); we defined this phenotype as "pla- 
code defect." In 30% (16/52) of cases we observed bundles of 
OE-type and VNO-type axons either overshooting past the OB or 
taking a misguided route (arrows in Figure 4). We also observed 
lack or impairment of connection with the OB, as indicated by 
the absence of typical glomerulus structures or their disorganized 
position at the OB (asterisks in Figure 4). We collectively defined 
these phenotypes as "connectivity and glomeruli defect." None of 
these phenotypes were seen in control embryos. 

z-dlx5a is the fish ortholog of mammalian Dlx5, in fact the 
embryonic expression territory is similar (86), and its knock- 
down causes craniofacial and neuronal phenotypes resembling 
the D/x5 _/ ~ phenotype in mice (87, 88). We depleted z-dlx5a in 
zebrafish embryos using a combination of two MO, and exam- 
ined the organization of olfactory axons. Following MO injection, 
72 hpf we recovered about 50% of CFP+ embryos (95% of the 
control injected) and about 72% of Venus+ embryos (78% of 
the controls). In 45% of cases (of 80 examined) we observed 
OPL defects, while in 54% of cases we observed OE-type and 
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Table 1 | Best Dlx5 target gene, selected combining the profiling results with PWM-based site prediction and embryonic expression. 



Gene title Gene symbol log2.FC. Dlx site Express Score Notes 





Leucine-rich repeat and Ig 
domain containing 2 


Lingo2 


-1.3995539 


+ 


N Ep 


5 


Structure similar to other Receptor Tyrosine 
Kinases, such asTrk. Associated to higher risk of 
tremor and Parkinson. Lingol is a component of 
the NOGO-66 receptor and may play a role in 
neurite outgrowth and oligodendrocyte 
differentiation 


Leucine-rich repeat LGI 
family, member 1 


Lgi1 


-1.355411 


+ 


N Ep 


5 


Secreted molecule of the SLIT family, promotes 
formation of stress fibers. Inhibits cell movement 
and invasion. Enhances growth of neuronal 
processes on myelin-based substrates. Its 
receptor forms complexes with Adam22 


Leucine-rich repeat protein 1, 
neuronal 


Lrrnl 


-1.032383 


+ 


N Ep 


5 


Transmembrane protein of unclear function. 
Regulates neurite growth 


Ig superfamily containing 
leucine-rich repeat 2 


Islr2 


-0.9967503 


+ 


Not/migr cell 


4 


Also known as Linx, could be a ReceptorTyrosine 
Kinase evolutionarily related toTrk receptor. 
Modulates axon extension and guidance 


ST8 ot-/V-acetyl-neuraminide 
a-2,8-sialyltransferase VI 


St8siaVI 


-1.3472121 




N Ep 


4 


Sialo-transferase expressed by neurons, essential 
for surface functions during neurite growth and 
neuronal migration 


A kinase anchor protein 6 


Akap6 


-1.5186358 


+ 


N Ep 


5 


Protein Kinase A-anchoring proteins. Serves as 
scaffold to bring together PKA and PDE and 
coordinate the timing and intracellular localization 
of cAMP signaling. Also binds to- and 
modulates-signaling through ERK, MAPK, and 
PP2A 


Dual adaptor for 
phosphotyrosine and 
3-phosphoinositides 1 


Dappl 


-1.2094534 




N Ep 


4 


Signaling adapter molecule, coordinates timing 
and location of signaling by PIP3 and PIP2 with 
that of ERK. Also binds F-actin and Rac 


RIKEN cDNA 9330120H11 
gene 


93301 20HllRik 


-1.1589186 




N Ep 


4 


Also known as HOMER 2, present at 
post-synaptic density, involved in receptor 
clustering, trafficking, and in calcium homeostasis 




Synaptosomal-associated 
protein 25 


Snap25 


-1.3481758 


+ 


N Ep 


5 


Controls membrane trafficking and fusion at the 
growth cone and at the synapse. Implicated in 
neuroblast migration and neuritogenesis during 
development. Forms complex with p140CAP 
which also binds to p130 CAS 


y-Aminobutyric acid (GABA) A 
receptor, subunit p 2 


Gabrb2 


1.0063619 


+ 


N Ep 


5 


Receptor subunit for GABA. GABA-b receptors 
mediate signals inhibitory for olfactory axon 
elongation 


Receptor transporter protein 1 


Rtpl 


-1.6548021 




N Ep 


4 


Chaperon, required for the efficient translocation 
of OR molecules to the membrane. Interacts with 
the OR and with Homer 


RIKEN cDNA 9330120H11 
gene 


93301 20HllRik 


-1.1589186 




N Ep 


4 


Also known as Homer2, present at post-synaptic 
density, involved in regulation of calcium fluxes 
and homeostasis 



(Continued) 
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Table 1 | Continued 



Gene title 



(D) AXON-GLIA INTERACTION PROTEINS 

Fatty acid binding protein 7, Fabp7 -1.9620307 
brain 



Gene symbol log2.FC. Dlx site Express 

N Ep-Gliale 



Score Notes 



Ermin, ERM-like protein 



Ermn 



-1.5033487 



N Ep/Sust 



Ganglioside-induced Gdapl 
differentiation-associated- 
protein 1 

UDP Galactosyltransferase 8A Ugt8a 



-1.1935825 + N Ep/Res 



-1.1139671 + 



Known as BLBP in human. Controls surface 
functions that are required for axon-Schwann cell 
interaction. May be involved in peripheral axon 
elongation and regeneration 

Also known as Juxtanodin. Expressed in 
sustentacular cells, binds to F-actin and stabilizes 
the actin cytoskeleton. In the CNS promotes 
myelination 

Involved in the Charcot-Marie tooth disease, in 
particular those forms with axonal deficits. 
Cellular function unclear 

Important for the biosynthesis of galacto-lipids 
and in myelin formation 



Cyclic nucleotide gated Cnga2 
channel a 2 

Visinin-like 1 Vsnll 



-1.129349 + N Ep 



-1.3860936 + N Ep 



Regulate axon extension and glomerular 
formation. KO mice have behavioral defects 
possibly linked to olfactory functions 

Also known as GP2. Calcium-regulated guanylate 
cyclase transduction system. Play a role in 
adaptation. Inhibits the formation of cAMR May 
affect dendrite and growth cone arborization 



Genes are sub-divided in five general categories (A-E). 



VNO-type axons targeting abnormal regions of the head near the 
OBs, often overshooting past the OB (arrows in Figure 4). We 
also observed impaired axon-OB connections, as judged by the 
absence of glomeruli-like bundles or their disorganized position 
(asterisks in Figure 4). None of these phenotypes were seen in 
non-injected or control MO-injected embryos. Thus, the deple- 
tion of z-dlx5a causes defects that recapitulates key aspects of the 
Dlx5~'~ phenotype (14, 16, 47). 

Next we focused on the putative Dlx5 targets Lrrnl, Lingo2, Islrl, 
St8siaVI, and Homer2, whose embryonic expression in the brain 
and olfactory system is reported in Figure S5 in Supplementary 
Material. We depleted z-lrrnl in 1-cell zygotes by MO injection. 
Of the injected embryos, only approximately 50% were recovered 
and positive for OMP:CFP (vs. 85% in the control injected), and in 
a majority of these (75% of a total of 62 examined) we observed a 
reduction of the CFP+ signal intensity. On the contrary, we recov- 
ered a not significantly different percentage of Venus+ embryos 
(71 vs. 78% in the control injected) and these occasionally (20%) 
showed a reduced YFP fluorescent signal (Figure 5). Twenty per- 
cent of z-lrrnl MO-injected embryos displayed placode defects, 
consisting in a reduced size, altered shape, and mispositioned neu- 
ron. Half of the z-lrrnl MO-injected embryos displayed an altered 
pattern of olfactory axon fasciculation and extension, with axons 
overshooting or taking an ectopic route (arrows) and reduced or 
absent glomeruli (asterisks). Thus, the depletion of z-lrrnl results 
in a delayed differentiation of the OMP+ type olfactory neurons 



and altered olfactory axons trajectory and connectivity, hallmarks 
of the phenotypes observed in Dlx5~ l_ mice and in z-dlx5a fish 
morphants. 

Next we depleted z-lingo2 in reporter zebrafish embryos. Injec- 
tion of the anti-z-lingo2 MO in 1-cell embryos caused minor OP 
defects, consisting in altered organization and shape, while axon 
trajectory and glomeruli formation appeared normal (Figure 6). 
Next we depleted z-st8SiaVI in zebrafish embryos. Injection of 
the anti- z-st8siaVI MO in 1-cell embryos resulted in a pheno- 
type affecting axon extension, trajectory, and glomeruli formation 
(Figure 6). Next we depleted Homer2 in the reporter fish embryos. 
Injection of the anti-z-homer2 MO in 1-cell embryos resulted in 
defects of OP organization and axonal targeting, plus also resulted 
in a reduced expression of Trpc2, seen as reduced YFP fluorescent 
signal (Figure 6). This last result might indicated that z-homer2 is 
involved in the differentiation of the VNO-type neurons, and its 
depletion may delay this process. Finally, the depletion of z-islrl 
yielded no appreciable phenotype affecting the olfactory pathway 
(data not shown). This gene is not expressed in the embryonic OE 
(Figure S5 in Supplementary Material). 

TESTING Dlx5, Dlx5 TARGETS, AND KS GENES IN ZEBRAFISH STRAINS: 
THE GnRH NEURONS 

To determine whether some of the DEGs that emerged from 
transcription profiling of Dlx5 mutants had some function of 
GnRH neuronal migration and neurite organization, we used the 
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FIGURE 4 | Depletion of endogenous z-dlx5a and z-fgfr1a in zebrafish 
embryos, to image the olfactory axons. (A) Micrographs of Trpc2:Venus 
(YFR yellow fluorescence) and OMP:CFP (cyan fluorescence) fish embryos 
injected with a control MO (top panels), injected with anti-z-d/x5a MO (middle 
panels) or injected with anti-z-fgfr7a MO (bottom panel). White arrows and 
lines indicate the normal axonal pathway in control embryos. Red asterisks 
indicate absence of glomeruli. Red arrows indicate altered axonal trajectories. 
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(B) Whole-mount bright field micrographs of injected embryo, showing an 
overall normal embryonic morphology and growth rate in the injected 
embryos, compared to control injected ones. (C) Proportions of embryos 
showing either placode defects (OPL disorganization, altered neuron 
distribution), or connectivity/glomeruli defects (altered axon trajectory, altered 
fasciculation, reduced or absent glomeruli), or both, upon injection of control 
(open bars), ani\-z-dlx5a (gray bars), or anti-z-fgfrl a (solid black bars) MOs. 
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FIGURE 5 | Depletion of endogenous z-lrrn 1 in zebrafish embryos, to 
image the olfactory axons. (A) Micrographs of Trpc2:Venus (YFR yellow 
fluorescence) and OMP:CFP (cyan fluorescence) zebrafish embryos 
injected with control (top panels) or with anii-z-lrrnl (bottom panels) MOs. 
The control MO did not cause any significant alteration. White arrows 
indicate the normal axonal pathway and glomeruli in the control embryos. 



Red asterisks indicate absence of glomeruli. Red arrows indicate altered 
axonal trajectories. (B) Whole-mount bright field micrographs of injected 
embryo, showing normal embryonic morphology and growth rate. 
(C) Proportions of embryos showing either OPL disorganization, or 
olfactory axon mistargeting, or both (last bars) upon injection of control 
(open bars) or ani\-z-lrrnl (gray bars) MOs. 



GnRH3:GFP transgenic zebrafish strain, previously reported (62- 
64). In these animals the GFP reporter is expressed under the tran- 
scription control of a fragment of the fish GnRH3 promoter. The 
G/iRFf3-GFP+ neurons have been widely characterized, and they 
consist in a population of terminal nerve associated GnRH+ neu- 
rons, thought to represent the mammalian hypothalamic neurons 
with olfactory origin (27, 62-64, 89) (Figures 7A,B). 

We depleted z-dlx5a, z-fgfrla/b, and z-lrrnl in the GnRH3:GFP 
1-cell zygotes, and examined the effect on the number, position, 
neurite organization, and commissure formation of the GFP3+ 
neurons associated to the terminal nerves. The depletion of z-dlx5a 
resulted in a reduced number of GFP+ neurons in 30% of cases, 
and in 70% of cases clearly appeared mispositioned (40 morphants 
examined) (Figures 7C,D; quantifications in 7G). However, the 
depletion of z-dlx5a did not affect the ability of GFP+ axons to 
cross the midline at the anterior commissure. Thus, a reduction 
of z-dlx5a in the fish model recapitulates (some of) the GnRH 
phenotype observed in the mouse model (16, 47). 

The depletion of z-fgfrla/b resulted in a reduced number 
of GFP+ neurons in 80% of cases (a total of 40 morphants 
examined), and in 22% of cases these neurons were clearly 



mispositioned, and had shorter neurites (Figure 7E). In 35% of 
cases the GnRH+ neurites failed to properly cross the midline in 
the anterior commissure. This phenotype recapitulates that seen 
upon depletion of z-kalla/b (45, 46), thus we conclude that, based 
on the results of two well-established KS/nCHH genes and one 
KS-causing gene in the mouse, the use of MO in the GnRH3:GFP 
strain is a valid approach to examine the KS phenotype in vivo, and 
assures that future analyses on this subject will be informative. 

The depletion of z-lrrnl in the GnRH3:GFP fish embryo caused 
a reduction in the number of GFP+ neurons in 45% of cases, 
and in 65% of cases caused their misposition along the terminal 
nerve (40 morphants examined) (Figure 7F). We also observed 
a reduction of their neurite length, but little of no defect of the 
anterior commissure. Thus also one Dlx5 target is involved in the 
organization and the maturation of the olfactory/GnRH system. 

BI0INF0RMATIC PREDICTION/PRIORITIZATION OF NEW KALLMANN 
DISEASE GENES 

A large set of genes has been found mutated, alone or in com- 
bination, in KS/nCHH patients, by classical mutation search 
approaches. However, with the exception of KALI and FGFR1, 
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FIGURE 6 | Olfactory and VNO axons, upon depletion of endogenous 
z-St8siaVI, z-lingo2, and z-homer2 in zebrafish embryos. (A) Micrographs 
of Trpc2:Venus (yellow fluorescence) and OMP:CFP (cyan fluorescence) 
embryos injected with control MO (top panels), or injected with 
ani\-z-st8siaVI, ani\-z-homer2, ant\-z-dapp1 , and ani\-z-lingo2 MOs, as 
indicated on top of each image. The control MO did not cause significant 
alterations. Arrows indicate altered axonal trajectory, asterisks indicate 



absence of glomeruli or altered OPL organization. Asterisks indicate the 
regions of reduced fluorescence intensity. (B) Whole-mount bright field 
micrographs of injected embryo, showing normal morphology and growth 
rate. (C) Proportions of embryos showing either placode defects, connectivity 
and glomeruli defects, or both, upon injection of the MOs indicated above 
(colored bars), compared to control MO (open bars). Asterisks indicate 
statistical significance. 
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FIGURE 7 | The GnRH3 neurons upon depletion of endogenous 
z-dlx5a, z-fgfr1a, and z-lrrn1 in zebrafish embryos. (A) Scheme 
showing the positions of the GnRH3:GFP+ neurons (green cells), 
relative to the OPL, the OB, and the olfactory nerves (yellow and blue), 
in a frontal view. The anterior commissure is shown at the basis of the 
OB. On top, a scheme illustrating the view plane (frontal) used in (A). 
(B) Scheme illustrating the view plane (ventral) used for the fluorescent 
images in (C-F). (C-F) Micrographs of GnRH3:GFP zebrafish embryos, 
at 60 hpf, injected at the 1-cell stage with either a control MO (C), with 



each of them is mutated in a small fraction of the patients, and 
together account for no more than 40% of KS/nCHH cases. Five 
novel genes, functionally linked to FGF8, have been recently iden- 
tified using predictive bioinformatics followed by mutation search 
in patients' DNAs (7). With the exception of some genes evidently 
linked (Prok2 and ProkR2; FGF8 and FGFR1; GnRH and GnRH- 
R) the remaining genes appear to be unrelated, or distantly related 
on a functional basis. We reasoned that relationships might exist 
between the KS-disease genes that are not obvious, or not easily 
detected, or that genes may have pleiotropic functions, not know 
as yet. Tools have been developed that search for such relationships 
in databases or newly generated data, and can be used to propose 
candidate disease genes (90). 

HUMAN NETWORK 

We compiled a list of genes known to cause KS, or KS and 
nCHH, excluding those causing only nCHH; the list included 
FGFR1, FGF8, KALI, PROK-2, PROKR2, CHD7, GnRH, GnRH- 
R, HS6ST1, TAC3, TACR3, SOX10 e SEMA3a. We also included 
FLRT3, IL17RD, FGF17, SPRY4, DUSP6, members of the "FGF8 
synexpression" group (7) and named all these "human reference 
genes." First we searched for the presence of the reference genes 
among the DEGs from the D\x5~ l ~ OE vs. WT, however none 
of them was found. Likewise, we searched for the presence of 
these genes among the DEGs from the time course of the nor- 
mal OE and VNO. With the exception of GnRH, none of the 



ar\i\-z-dlx5a MO (D), with ar\i\-z-fgfrla MO (E), or with ar\t\-z-lrrn1 MO 
(F). Insets on the lower right of each panel is a low-magnification 
merged micrograph (bright field and GFP fluorescence) of the higher 
magnification one. Red asterisks indicate reduced number of cells, red 
arrowheads indicate scattering and delayed cell migration, white 
arrowheads indicate the anterior commissure, white asterisks indicate 
absence of anterior commissure. (G) Quantification of the observed 
phenotypes, as percent over the total number of GFP+ embryos 
examined with each MO. 



other genes was found. Next, we positioned the "human refer- 
ence genes" within the global conserved co- expression network, 
using TS-CoExp (with the exception of KALl/anosminl that lacks 
a mouse ortholog and for which the conservation criterium cannot 
be applied) (Figure 8A). 

The network representation does not consent per se to derive 
relevant information. Instead, from the data we extracted those 
genes connected with at least six (N = 2), at least five (N = 3), at 
least four (N = 10), at least three (N = 45), at least two (N = 317), 
and at least one (N = 1977) reference genes. We then categorized 
these genes by G.O. and detected an enrichment of the following 
G.O. categories: phosphoproteins, kinase/transmembrane recep- 
tors, cell adhesion, cell junctions, regulation of cytoskeleton, cell 
migration/motility, neuronal projection. Among the most con- 
nected ones we did not find any gene causing KS in mice, but we 
identify TRIM2, GATAD2A, SNRPN, and CDH2. Being expressed 
in the embryonic OE (Figure S6 in Supplementary Material), these 
represent most interesting genes. 

Next, we used the disease gene prediction tool of TS-CoExp 
to identify novel candidate KS genes: the "human input genes" 
were taken as reference to independently prioritize the follow- 
ing DEG lists: (a) Dlx5~ l - OE vs. WT at E12; (b) EPI OE 14 vs. 
OPL Ell (WT); (c) EPI VNO 14 vs. OPL Ell (WT); (d) MES 
OE14 vs. OPL Ell (WT); (e) MES VNO 14 vs. OPL Ell (WT). 
From the DEG list (a) we found nine genes significantly associ- 
ated with the KS phenotype, three of which (RGS5, F2RL1, and 
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DPF3) are expressed in the embryonic OE, while two (GATA3 and 
ADAMTS5) are expressed in the olfactory mesenchyme (Figure S6 
in Supplementary Material). 

From the DEG lists (b) and (c) we found 2 1 and 73 genes respec- 
tively, significantly associated with the KS phenotype, 19 of which 
are present in both lists, and the majority of these are expressed in 
the nasal mesenchymal. Notably, the search predicted two genes 
known to cause KS in the mouse, namely Ebf2 and Nrpl, con- 
firming that our analysis is in principle correct. From the DEG 
lists (d) and (e) we found 47 and 189 genes respectively, signifi- 
cantly associated with the KS phenotype, 27 of which are present 
in both lists. Among these, 50% show expression in the embry- 
onic OE (ACAN, AKAP6, ATF5, KRT18, MYTH, NDRG1, NRXN1, 
SYT1, and TPD52) and 50% in the olfactory- associated mes- 
enchyme (ANXA1, DCN y FCGRT, PAPS2, PTRF, RUNX1, SlOOb, 
and TGM2). Notably, the search predicted two genes known to 
cause KS in the mouse, namely Ebf2 and Nrpl. Furthermore, 
we found genes such as AKAP6, LING02, LGI1, and LGI2 that 
were found among the Dlx5 targets in the OE, and SEMA3C and 



MET, known to play a role in axon guidance and cell migration, 
respectively. 

MOUSE NETWORK 

We applied the approach previously used to those genes causing 
a KS-like phenotype in mice; the list included Dlx5, Emx2, Klf7, 
Fezfl, Sixl, Prok2, ProkR2, Lhx2, Shepl, Ebf2, Nrpl, and Sema3a. 
None of them was found in the DEG lists from the Dlx5 _l ~ OE 
vs. wild-type, and none was found in the DEG list from the time 
course of the normal OE and VNO, with the exception of Lhx2 
and GnRH. We positioned the mouse reference genes within the 
conserved co- expression network (with the exception of Fezfl for 
which no result was obtained) (Figure 8B), and extracted lists of 
genes connected with at least five (N = 1), at least four (N = 3), 
at least three (N = 33), at least two (N = 261), and at least one 
(N = 1850) of them. These genes were then categorized by G.O. 
and we detected these over- represented terms: phosphoproteins, 
cell adhesion/cell junctions, neuronal projection, cell motility, 
cytoskeleton regulation, transcription regulation. 




FIGURE 8 | Continued 
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FIGURE 8 | Disease gene networks for KS. (A) Position of the human these analyses are in the text. For simplicity, only the genes connected with 

KS-causing genes within the global conserved co-expression network, as at least three input genes are shown; the genes connected with "at least 

computed with theTS-CoExp algorithm. (B) Position of the genes causing a one" or "at least two" input genes are available upon request. Green circles 

KS-like phenotype in mice, within the global conserved co-expression represent the input genes, pink circles represent the connected genes, lines 

network. The lists of the input (human and mouse) "disease" genes used for represent statistically significant co-regulations. 
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Next, for each of the four DEG list (from the profiling results, 
see above), we used the disease gene prediction tool in TS-CoExp 
to identify those genes most likely to be involved in KS, using the 
mouse reference genes as input. From the (a) list we found 13 
genes, two of which (Scn3B and Sv2B) are expressed in the OE, 
while two (Adamts5 and Wnt5a) are expressed in the olfactory 
mesenchyme (Figure S6 in Supplementary Material). From the (b) 
and (c) lists we found 19 and 65 genes respectively, 17 of which are 
present in both lists. Most are expressed in the nasal mesenchyme. 
Contrary to the human network, the mouse network did not pre- 
dict any human KS gene. From the (d) and (e) lists we found 65 
and 41 genes respectively, 17 of which are present in both lists and 
most of these are expressed in the nasal mesenchyme. Also in this 
case, we could not predicted any human KS genes. We found two 
genes: Dcx and Intra, the first is relevant for migration of imma- 
ture neurons, the second codes for a leucine-rich repeat protein 
similar to Lrrnl. 

Finally, considering both the human and the mouse reference 
genes, five genes were found in at least three lists of genes associated 
with KS, namely Den (Decorin), FGF7, Aspn, Ptfr, and Ntrk2. Three 
of these, Ntrk2, Den, and Ptfr are clearly expressed in the juxta- 
OE and VNO mesenchyme, Aspn is ubiquitously expressed and 
the expression of FGF7 is unclear. FGF7 codes for a growth factor 
related to FGF22, the literature reports indications that it functions 
as pre-synaptic organizing molecule during hippocampal develop- 
ment (91, 92) and is needed for migration of enteric neurons (93). 
Its function in the developing OE is unknown. In conclusion, the 
most promising predicted KS genes are mesenchymally expressed. 
This is not surprising, considering that in all our profiling results, 
the prevalent categories are cell-cell and cell-matrix interactions, 
remodeling, signaling, etc. . . 

Note: all gene lists, categorizations or tables not included in the 
manuscript or as supplementary material are available on request. 
All data are deposited at GEO repository, N° GSE52800. 

DISCUSSION 

Kallmann syndrome and nCHH are developmental/pediatric con- 
ditions phenotypically well characterized, however less well under- 
stood molecularly. Despite the number of genes found mutated in 
KS/nCHH patients, the majority of them still await a molecu- 
lar definition. Thus, there is a strong basis to predict that many 
additional disease loci remain to be identified. Furthermore, the 
mutations found in KS patients, once thought to act alone, are now 
recognized as cooperating mutations, and the prevalent notion 
states that most KS/nCHH cases should be a bi-genic or oligo- 
genic disease (21, 22). This raises hopes that a more exhaustive 
knowledge of cooperating genes and mutations, should consent a 
better prognosis and possibly personalized therapies. 

Methods and algorithms have been proposed to iden- 
tify/prioritize novel disease genes, based on (meta) -analyses of 
specific profiling data, co-expression networks, genome locations, 
functional categorizations, protein -protein interactions, etc. (90, 
94). These methods have several advantages over direct whole 
exome sequencing of large panels of DNAs from KS/nCHH 
patients (95-97). In this study, on one side we have uncov- 
ered functional classes, possible networks, and individual genes 
involved in the olfactory/GnRH developmental, and validated 



some of them in the zebrafish model recapitulating the KS pheno- 
type. On the other side we positioned known human genes causing 
KS and mouse models with a KS-like phenotype in gene-co- 
expression networks, in order to identify genes potentially relevant 
for the process and candidate KS-disease genes. 

Embryonic development of the olfactory connection and the 
migration of immature GnRH neurons are anatomically and func- 
tionally linked. Since olfactory detection is a primary sensory 
system in most vertebrates, and sexual maturation/reproduction 
is essential for the species, it is not surprising the developmental 
process is highly conserved and is controlled by multiple - partially 
redundant - networks of molecular regulations. The high degree 
of conservation among vertebrates justifies the use of the zebrafish 
embryo for in vivo testing (98): not only its general anatomy has 
not greatly changed, but also the migration of GnRH neurons 
along the terminal nerve, in association with the VNO axons has 
been overall maintained (83). 

We have generated profiling data, comparing the mouse normal 
OE and VNO at three developmental time-points, and comparing 
the normal vs. a mutant model characterized by a KS-like pheno- 
type. The data have been used to identify novel gene categories 
involved in the development of the olfactory system, to identify 
Dlx5 target genes in the OE, and to intersect this wealth of infor- 
mation with data from other sources. As a further step, it might be 
useful to generate profiles from other models of KS in the mouse 
(i.e., Prok2, Fezfl y etc.) and intersect the results searching for com- 
mon patterns and co- regulations. We have attempted this, however 
with little success, most likely because we specifically profiled the 
olfactory epithelia, while data form the Emx2 and Klf7 models 
were generated from the OB. Likewise, it would be useful to inter- 
sect our profiles from the Dlx5 model with datasets from freshly 
dissociated embryonic GnRH neurons. 

Categories that emerge from the "time course" profiles strongly 
implicate extracellular matrix remodeling, cell adhesion, and cell- 
cell signaling molecules. This is true both for the OE and VNO 
development, that after all appear more "similar" than "different." 
The profiles of the "pathologic" condition, i.e., the Dlx5 knock-out 
model, indentified a number specific molecules in the categories 
of membrane receptors/adhesion molecules, axon-glia interaction 
molecules, but nothing specifically related to "axon elongation." 
This suggests that the cell- autonomous properties of the olfac- 
tory/VNO axons to establish connections, provided (directly or 
indirectly) by the transcription factor Dlx5, reflects cytoskeletal 
properties and cell surface events, mediated by receptors, scaffold 
proteins and cell adhesion (see below). 

Since profile data may easily lead to false positives, functional 
validations are mandatory; we show that the zebrafish embryo 
can be effectively used either to examine the trajectory and of the 
olfactory axons, or the status of the GnRH3 neurons. We have 
functionally tested five genes for olfactory axons, and three genes 
for GnRH3 neurons, and the results clearly indicate that the chosen 
genes do affect axonal trajectory and GnRH3 migration. Previous 
works have shown that the depletion of z-kalla/kallb in the fish 
embryo also causes KS-like phenotypes (45, 46), thus the use of 
the reporter zebrafish strains we have adopted appears to be a valid 
approach in which to examine new KS-causing genes in human, 
or genes causing a KS-like phenotype when mutated in mice. 
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Work is continuing in this direction. The following interesting 
genes/categories emerge from the profile data. 

THE LEUCINE-RICH REPEAT PROTEINS 

We identified three leucine-rich repeat transmembrane protein 
genes among the Dlx5~ /_ targets, namely Lrrnl, Lingo2, and Lgil. 
We functionally tested two of these using zebrafish embryos, and 
the results clearly show that these proteins participate in the 
development of the olfactory pathway. Lrrnl was also tested in 
GnRH3:GFP fish embryos, and the results show that it is required 
for correct GnRH neuron migration. Furthermore, Lingo 1, Lingo2, 
and Lgi were prioritized with the human network, and Lrrtm 
(another member of this family) emerged from the mouse 
network. 

Lrrnl is a glycosylated single-pass transmembrane protein 
with 12 external leucine-rich repeats, a fibronectin domain, an 
immunoglobulin domain and short intracellular tails capable of 
mediating protein-protein interaction. Lrrnl is closely related 
to drosophila tartan/capricious (trn/caps) proteins. Differential 
expression of trn/caps promotes an affinity difference and bound- 
ary formation between adjacent compartments in a number of 
contexts. The regulated embryonic expression and cellular location 
of these proteins suggest important roles during mouse develop- 
ment in the control of cell adhesion, movement, or signaling (99). 
Indeed, Lrrnl has been identified as a positive and negative regula- 
tor of neurite growth (100). Lrrnl appears to be a key regulator of 
the process of generating distinct cells at the midbrain-hindbrain 
boundary of the brain. In the chick embryo Lrrnl is dynami- 
cally expressed, the timing of its down -regulation correlates closely 
with the activation of signaling molecule expression at bound- 
ary regions. Cells over- expressing Lrrnl violate the boundary and 
this result in a loss of cell restriction at the midbrain-hindbrain 
boundary (101). Lrrnl may regulate the subcellular localization of 
specific components of signaling or cell-cell recognition pathways 
in neuroepithelial cells (102). 

Lingo2 is an exclusively neuronal transmembrane pro- 
tein (103), containing 12 extracellular leucine-rich repeats, an 
immunoglobulin C2 domain and a short intracellular tail, and 
with a predicted structure similar to the Trk Receptor Tyrosine 
Kinases. In human Lingo2 been linked both to essential tremor 
and to Parkinsons disease (104). Interestingly, the combination of 
leucine-rich repeat and immunoglobulin-like domains is found in 
the domain architecture of the Trk neurotrophin receptor protein. 
In the mouse embryo Lingo2 is expressed in a the olfactory neurep- 
ithelium and in various areas of the adult brain (99). Lingo 1, 
another neuron-specific member of the same family, has been 
shown to be a component of the Nogo66 receptor/p75 signal- 
ing complex (74). This ternary complex confers responsiveness to 
oligodendrocyte myelin glycoprotein, as measured by RhoA acti- 
vation. Such responsiveness is linked to the inhibition of axon 
regeneration of neurons in the adult brain, by myelin. Thus, 
Lingo proteins are likely to play a role in neurite outgrowth and 
oligodendrocyte differentiation. 

Lgil is a leucine-rich repeat molecule, found to be down- 
regulated in the absence of Dlx5. This is a secreted molecule 
of the SLIT family, promotes formation of stress fibers, inhibits 
cell movement and invasion, and enhances growth of neuronal 



processes on myelin-based substrates (75, 76). At the moment we 
have no functional data on the possible role of Lgil in olfactory 
development, yet should be explored. 

MULTI-ADAPTOR SCAFFOLD MOLECULES 

Among the Dlx5 targets we note the presence of a set of scaffold- 
adaptor proteins, including Akap6, Dappl (also known as BAM32), 
and Homer2. Akap6 belongs to a class of protein kinase A- 
anchoring proteins, serving as scaffolds to cluster PKA and PDE 
and to coordinate the timing/intracellular localization of cAMP 
signaling. Akap proteins also bind to- and modulate-signaling 
through ERK, MAPK, and PP2A (105, 106). The potential impor- 
tance of this class of molecules is suggested by the fact that Akap6 
(expressed in olfactory neurons) and Akap2 (expressed in the mes- 
enchyme) emerge as predicted/prioritized disease genes from the 
human network. Akap6 is absent in the zebrafish genome and 
could not be tested. 

Dappl codes for a signaling adapter molecule, much studied 
in B lymphocyte activation, in which it coordinates timing and 
location of signaling by PIP3 and PIP2 with that of ERK. Dappl 
also binds F-actin and Rac (107-109). Dappl is not apparently 
expressed in the embryonic OE, nevertheless when Dappl was 
depleted in the fish model a mild effect on axonal trajectory and 
OPL organization have been observed. It appears very likely that 
lipid signaling is involved in axonal trajectory and connectivity 
during olfactory development. 

Homer2 is a post-synaptic scaffold molecule, involved in recep- 
tor clustering and modulating their downstream signaling. How- 
ever, recently a role for Homer2 in tuning the activity of G- 
proteins coupled receptors (such as the ORs) by controlling cal- 
cium influxes has been demonstrated. We carried out functional 
experiments depleting Homer2 in zebrafish embryos: the results 
provide evidence for its involvement in olfactory axonal devel- 
opment. Considering the established importance of the OR for 
olfactory axon connectivity and guidance during embryonic devel- 
opment, much before their role in odor perception, an embryonic 
role of Homer2 can be envisioned, and our results with zebrafish 
clearly show this. 

The p130CAS - Shepl regulation 

Mouse embryos null for Shepl show retarded OE differentiation, 
lack of primary axonal connections with the OB and retention of 
GnRH neurons in the nasal mesenchyme (110). These defects are 
accompanied by a reduced phosphorylation of the multi- adapter 
scaffold molecule pl30CAS in the olfactory neurons and axons. 
Shepl promotes Src-dependent phosphorylation of the multi- 
adapter molecule pl30CAS, in vitro (111). These data implicate 
the phosphorylation of pl30CAS in the establishment of olfactory 
contacts and in GnRH neuron migration, in line with previous 
studies suggesting that pl30CAS is required for neurite outgrowth 
and axon guidance (112-114). pl30CAS belongs to a family of 
multi-adaptor and scaffold molecules that spatially and tempo- 
rally collect, integrate, and modulate signals coming from RTKs 
and adhesion receptors (115-117), undergoing changes in phos- 
phorylation and interacting with a large set of effectors proteins. 
In light of the phenotype of Shepl' 1- mice, the involvement of 
pl30CAS in olfactory development and GnRH neuron migration 
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is a likely possibility to be explored. Since pl30CAS null mice 
are embryonic lethal (118), this study will have to be pursued via 
conditional deletion of pl30CAS in the olfactory system. 

St8siaVI codes for a sialyl-transferase, expressed by olfactory 
neurons. The highly related St8siaII and St8siaIV proteins are 
required for polysialylation the N-CAM, confer to this neuronal 
surface molecule an ti- adhesive properties and thereby promote 
neurite elongation and cell migration (77-79). Thus a role for this 
protein in OE development is conceivable, and supported by the 
presented data in fish embryos. 

EphA3 codes for a receptor for the guidance molecules Eph- 
rinA3 and EphrinA5, which are expressed by VNO axons and have 
a preference for interacting with EphA expressing cells in the Acces- 
sory OB. Alterations of this pathway leads to abnormal topography, 
i.e., guidance defects, of the olfactory and VNO axons (119), indi- 
cating the EphrinA-EphA system is a positive guidance cue. Dlx5 
is co-expressed with EphrinA3 and EphrinA5 in the VNO, while 
EphA3 is expressed in the mesenchyme near the VNO (Figure S7 
in Supplementary Material). The link between Dlx5 and EphA 
signaling should be deeply explored. 

GENES EMERGING FROM THE BIOINFORMATIC ANALYSES 

A recent work has succeeded to use bioinformatics to priori- 
tize candidate KS genes, focusing on the FGF8 co- expression and 
functional network (7). Inspired by this work, we opted for an un- 
biased approach that simultaneously searches for links between 
genes apparently unrelated. Limiting our search to co- regulations, 
we strongly introduce the notion of conservation, reasoning that 
the olfactory/GnRH development is highly conserved within ver- 
tebrates. Indeed, in our work we have attempted to use also mouse 
KS-disease genes to run the search. The advantage is the pos- 
sibility to use all the disease genes, instead of focusing only on 
those logically related. An additional advantage of the present 
work derives from combining bioinformatic predictions, putative 
gene functions, phenotype descriptions, and information from 
the literature with "wet" profiling data specifically obtained from 
embryonic olfactory tissues. 

The "human" network was able to predict few mouse KS gene 
(Ebf2 and Nrp J), providing an evidence that the algorithm is effec- 
tive. The outcome, both in terms of individual genes and the G.O. 
classes, assures that the pipeline works. The addition of protein- 
protein interaction data (when made available) or other data to 
carry out meta-analyses will certainly refine the results. On the 
contrary, it appears that the "mouse" gene network is little infor- 
mative, i.e., less able to predict the human KS genes. This might 
be due to the fact that the definition of mouse input gene is based 
on accurate phenotypic analyses on the olfactory system, reported 
in the literature, that scientist don't routinely conduct (we may 
miss many other genes) or it is incomplete and does not examine 
olfactory axons but only hypothalamic GnRH neurons. 

The predicted/prioritized genes emerging from our analyses 
may represent a novel set of KS-causing genes, or genes that might 
contribute when co-mutated with others. While the use of modern 
sequencing approach (WES) on KS patients' DNAs is the straight- 
forward approach to define their role in the human disease, addi- 
tional filters may be needed to further prioritize these genes, i.e., 
testing their function on GnRH3+ neurons fish embryos. 



TRIM2 - tripartite motif containing 2, codes for an E3 ubiquitin- 
protein ligase that has been implicated in ubiquitination of neu- 
rofilament light chains. TRIM2 controls the dynamic of neu- 
ronal cytoskeleton, by which determines the specification of the 
choice of the axonal vs. dendritic projection in hippocampal 
neurons (120). 

CDH2 - cadherin-2, also known as N-cadherin, codes for a 
well known calcium-dependent neuronal cell adhesion molecule 
that contributes to the formation of neural circuits by regulating 
growth cone migration and synapse formation. In the mam- 
malian embryonic neocortex, radial migration is instructed by 
several signals that include homophilic interactions mediated by 
Cdh2 (121), and the fish embryo Cdh2 is involved in neurob- 
last migration within the hindbrain (122, 123). Chd2 function 
is required for guidance of afferent fibers of cranial sensory 
neurons (124) and regulates motor axon growth and branch- 
ing, in fish embryos (125). During olfactory development, Cdh2 
is expressed by receptor neurons and closely parallels expres- 
sion of y-catenin in neuronal axons (126), thus Cdh2 is posi- 
tioned to underlie the formation of olfactory primary olfactory 
connections. 

ADAMTS5 codes for a disintegrin-like and metallopeptidase 
extracellular protease, with thrombospondin-like motif. Adamts5 
plays a role in the specification and patterning of progenitor cells 
in the lateral and medial ganglionic eminences (127). The pro- 
teolytic cleavage of astrocyte-derived proteoglycan, exerted by 
Adamts5, loosens the matrix environment and promotes neu- 
rite outgrowth (128). Being predicted by both the human and 
the mouse disease-gene networks, Adamts5 appears to be a very 
interesting candidate. 

RGS5 - regulator ofG-protein signaling 5, codes for a protein that 
accelerates the inactivation of Ga-dependent signaling in various 
cells types. Down-regulation of RGS5 induces GPCR-mediated 
signaling pathways and promotes migration of vascular and cancer 
cells (129, 130). A role of this protein in promoting the migration 
of GnRH neurons is possible, although RGS5 null mice don't show 
any obvious phenotype (131). 

DPF3 - D4, zinc and double PHD fingers, family 3, codes for a 
component of the BAF chromatin remodeling protein, and acts a 
transcription co-activator in SWI/SNF complex- activation (132). 
DPF3 functions to activate transcription of the target genes Pitx2 
and Jmjdl c in association with the BAF complex, and binds histone 
H3 and H4 in an acetylation-dependent manner (133, 134). How 
this could be relevant for olfactory development, GnRH neuron 
migration and KS, is unclear. 

FGF7 has been proposed to act as a pre-synaptic organiz- 
ing molecule in the mammalian brain, and in particular dur- 
ing hippocampal development. Indeed FGF7-deficiency impairs 
inhibitory synapse formation, which results in mossy fiber sprout- 
ing and enhanced neurogenesis (91, 92). Neutralization of FGF7 
inhibits pre-synaptic differentiation of mossy fibers at contact with 
granule cells, and inactivation of FGFR2 has similar effects (92). 
In neurons, FGFs and cell adhesion molecules stimulate neurite 
outgrowth via activation of FGF receptors. A role for FGF7 for 
the migration of enteric neuroblasts has been suggested from 
analyses of CAMs and FGFs expression in Hirschsprung Disease 
patients (93). 
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CONCLUSION 

The molecular control over the ability of olfactory axons to con- 
tact the anterior forebrain, and/or the ability of GnRH neurons to 
efficiently migrate and home to the hypothalamus, entails numer- 
ous proteins of various functional classes, many of which appear 
to be directly and indirectly involved in matrix remodeling and 
signaling. Indeed, the data indicate that the navigation of OE and 
VNO axons is mostly governed by cell-cell and cell-matrix cues, 
rather than intrinsic properties of the axons. These include a set 
of scaffold molecules that, for their nature, are strong candidates 
for playing a key role in guiding axonal elongation-guidance and 
connectivity, as well as for GnRH neuron migration and hom- 
ing. These molecules will be of great interest for developmental 
biologists. 

Perturbations in the expression and sequence (mutations) of 
these molecules and in their associated gene networks may cause 
phenotypes similar to KS, a possibility that can be rapidly tested 
in zebrafish strains, and eventually in the mouse. Human geneti- 
cists should consider these molecules for mutation screens. This 
opens the possibility to test them in the mammalian model and to 
search for mutations in large collections of DNAs from KS/nCHH 
patients, hereditary, or sporadic, with the hope to find mutations, 
alone or in combination with mutations in known KS/nCHH 
genes. 

Finally we show the validity of approaches based on high- 
throughput data generation and predictive bioinformatics to iden- 
tify genes potentially relevant for specific developmental processes, 
and ultimately for disease. Indeed, we have uncover a set of mole- 
cules that might be candidate disease genes, to be tested in future 
mutation screens. 
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